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Abstract 

We use numerical simulations to study the phase behavior of a system of purely repulsive soft 
dumbbells as a function of size ratio of the two components and their relative degree of deforma- 
bility. We find a plethora of different phases which includes most of the mesophases observed in 
self-assembly of block copolymers, but also crystalline structures formed by asymmetric, hard bi- 
nary mixtures. Our results detail the phenomenological behavior of these systems when softness is 
introduced in terms of two different classes of inter-particle interactions: (a) the elastic Hertz po- 
tential, that has a finite energy cost for complete overlap of any two components, and (b) a generic 
power-law repulsion with tunable exponent. We discuss how simple geometric arguments can be 
used to account for the large structural variety observed in these systems, and detail the similarities 
and differences in the phase behavior for the two classes of potentials under consideration. 
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I. INTRODUCTION 



The problem of optimal packing of spheres in two and three dimensions is a subject 
that has inspired multiple generations of scientists. The first conjectures for equally sized 
spheres were put forward by Kepler at the beginning of the seventeenth century, and were 
finally proven in 1998 by Thomas Hales pp. Far from being a merely academic problem, 
packing and self-assembly of nanoparticles and mesoscopic components have become major 
challenges in nanotechnology. Early experiments on colloidal dispersions [2] have shown how 
hard spherical particles can spontaneously organize into ordered structures when compressed 
at sufficiently large volume fractions, thus providing a microscopic realization of the spherical 
packing problem. The nature of this fluid to solid transition is understood in terms of the 
configurational entropy of the system which, at hight densities, is maximized when the 
particles reside on the lattice sites of an FCC or an HCP crystal. The relative free energy of 
FCC and HCP structures for hard spheres was studied in detail by computer simulations [3] . 

Although the specific case of hard spherical particles is a very relevant one, as it represents 
probably the simplest nontrivial mathematical reference model, the number of ways in which 
such particles aggregate into ordered structures is very limited; even in the presence of 
additional isotropic attractions due for instance to dispersion forces as in the Lennard- Jones 
interaction or to depletion forces as in colloid-polymer mixtures. 

Recent advances in particle synthesis |IHTT] have significantly extended the structural 
landscape via the development of colloidal particles that are anisotropic both in shape and 
surface chemistry, thus providing an unlimited number of building blocks that can spon- 
taneously organize into an unprecedented variety of structures with potentially novel func- 
tional, mechanical, and optical properties. For this reason, the problem of self-assembly of 
nanoparticles is today under intense investigation. Although most of the work has focused 
on how specific interactions (typically attractive) may be tailored to drive self-assembly of 
nanocomponents, over the past few years it has been shown how a large number of struc- 
tures with substantially different packing symmetries than FCC, HCP or BCC may become 
accessible even to purely repulsive spherical particles, either by altering the shape of their 
mutual interactions or by relaxing the constraint of mutual impenetrability. The formation 
of mesophases by particles interacting via a hardcore plus a repulsive shoulder potential was 
analyzed in [12j, and the relative stability of A15, BCC and FCC crystal structures as a func- 
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tion of particle softness was discussed in [13] . More complex mesoparticles such as charged 
or neutral star polymers, dendrimers or microgels are usually dealt with by developing a 
coarse-grained representation of the components which are reduced to objects with a simple 
shape interacting via ad-hoc effective pair potentials. The study of the physical properties of 
systems adopting these exotic intercations, some of which allowing for even complete overlap 
among the components, has been the subject of several publications [TH - 122] . These studies 
unveiled the existance of an unprecedented richness in phenomenological behavior including 
the possibility for some classes of interactions to form polymorphic cluster phases (23] or 
reentrant melting transitions followed by multiple transitions between ordered crystalline 
phases [2H - T26] . and set the stage for the work presented in this paper. For a recent review 
on the subject we refer the reader to reference [27] . 

In our recent work [26J, we have shown that a large number of solid phases becomes 
accessible upon increasing the density of a system of spherical particles interacting via a 
Hertz potential. Among the many phases we found crystals with cubic, trigonal, tetragonal, 
and hexagonal symmetries. In this paper we go one step further and include an extra degree 
of complexity to the problem: particle geometry. Specifically, we explore the relative role 
of softness and geometry in a system of asymmetric, repulsive dumbbells. We present the 
system phase behavior as a function of external pressure and discuss the interplay between 
geometry and softness as a means of understanding the formation of the large number of 
ordered phases found in the present work. We have chosen the asymmetric dumbbellar model 
because it allows us to efficiently and systematically incorporate both geometrical asymmetry 
and particle softness into the single components. In order to keep our results as general as 
possible we decided to describe the interparticle effective interactions using two common 
functional forms: (1) the Hertz potential, which describes the energy cost associated with 
elastic deformation of an amount h between two spheres, and has the form V(h) oc Yh 5 ^ 2 ( 
Y is the Young's Modulus) [2S], and (2) a generic inverse power-law V(r) oc r _n cut off at 
the particle diameter. 

In the first case we tune the softness by changing the Young's modulus Y, whereas in the 
second case this is done by changing the exponent n. 

When both tuning parameters tend to infinity, the two models become identical and the 
hard-repulsive limit is recovered. It is important to stress that unlike the inverse power 
law interaction, the Hertz potential is bound, i.e. there is a finite energy cost to completely 
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overlap two components. In this respect the two potentials can be considered as belonging to 
different classes of interactions, and we therefore discuss them independently. The schematic 



We begin by considering the elastic pair potentials. The dumbbells are modeled as 
two soft spheres of diameter a± and 02, connected with a rigid bond of equilibrium length 
d = 1/2(0"! + cr 2 ). The softness of the spheres is modeled via the Hertz potential, that 
describes the change in the elastic energy of two deformable spheres when subjected to axial 
compression [28, 29J. The potential has the following form: 



The indices i,j G 1,2 indicate the identity of the spheres, is the distance between the 
centers of any two spheres that belong to different dumbbells, and Yi is the Young's modulus 
of the sphere. In the equation above the Poisson ratio is taken to be 1/2. When YiYj/(Yi + 
Yj) — > 00 the interaction between the two spheres approaches the hard-core limit. When 
YjYj/ (Yi + Yj) —7- the spheres can completely overlap. The important property of the Hertz 
potential is that it is bounded and remains finite at = 0. 

We used the LAMMPS molecular dynamics package [30] with a Nose/Hoover thermostat 
in the NPT ensemble to study the phase behavior of a system of 512 dumbbells for a wide 
range of system densities, for different values of Y% and Y2 and several values of size ratio 
s = 02/01. The simulation box was set to be a cube with periodic boundaries for most of 
our simulations, no difference in the results was found when decoupling the box lengths in 
each of the three Cartesian coordinates. The system initial configurations were prepared by 
performing NPT simulations in the gas phase at a very low pressure p. Once the system 
is equilibrated, we slowly ramped the pressure to the desired final value pj starting form 
Po in small increments Ap. Each subsequent simulation performed at the constant pressure 
Pi started form the thermalized configuration at the pressure pj_i = Pi — Ap. This slow 
pressure annealing procedure should allow the system to fully equilibrate at each pressure. 
For each Ap, simulations were run for a minimum of 2 ■ 10 6 steps with the reduced time-step 
size St = 0.008to ( where to is the unit time expressed in standard MD units). 

Our system has four independent parameters: namely s, Yi, Y 2 and system pressure p, 



representation of our model for a soft asymmetric dumbbell is shown in Fig. [TJ while the 
pair potentials we use are illustrated in Fig. [2| 




(1) 
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making the phase diagram very complex and multidimensional. We find that a convenient 
way of presenting the results is by showing a sequence of slices in the {Yi, Y 2 } two-dimensional 
plane, each at constant s and p. To further reduce the size of the phase space we impose 
the following restrictions: Y\ > Y 2 and o\ < o 2 , i.e. the smaller component of radius 0\ is 
always the one that is harder to deform. Given our previous work [22], we also expect this 
constraint to result in the most interesting structure formation. 

Finally, all pressures have been rescaled to properly account for the different sizes of the 
dumbbellar components, p* = p(s 3 + Since there are no attractive forces in our 

system, the formation of these phases is mostly due to the minimization of the repulsion 
energy as a response to the applied external pressure. However, entropic contributions 
are not completely negligible [25] . especially in the case of dumbbell-crystallization at low 
pressures, as will be described below. 

Let us start by looking at the hard-limit behavior, i.e. when Y\ — Y 2 > 10 4 , depicted 
in Fig. [3] When the two components have equal diameter (and up to s = 1.05), the 
dumbbells, as expected, crystallize into an FCC crystal. We also find that the system can 
easily crystallize for size ratios in the range 2.4 < s < 3. In this case the dumbbells pack 
into a NaCl-type binary crystal. Interestingly, we find that for intermediate values of s 
(1.05 < s ^ 2.4) the system remains fluid, implying that the hard limit has a crystallization 
gap in the size ratio space. These results are a bit surprising as a few papers dealing with the 
equilibrium properties of binary AB-type hard-sphere mixtures and hard dumbbells [3"T143"3"] 
report the stability of a variety of crystalline structures across a comparable range of size 
ratios. One would expect some of these structures to also nucleate from the fluids phase in 
our system, especially the ones predicted for hard dumbbells [33]. We believe that this is 
only an apparent discrepancy, and can be resolved in terms of the kinetics of the system. 
In fact, it is well known that the probability of crystal nucleation of hard spheres becomes 
quickly reduced above a certain degree of polydispersity [5ij . It is therefore possible that the 
crystals reported in [33] are indeed the most stable ones. However, it is also possible that for 
1.05 < s ^ 2.4 they are kinetically inaccessible and for 2.4 < s < 3 the closest in free energy 
to the partent fluid phase is simply the NaCl-type. What is also very interesting is that as 
soon as the hard limit is abandoned, the crystal phase space increases in size and diversity: 
the upper limit of this crystallization gap is narrowed to s = 1.8, and other crystal structures, 
besides NaCl-type, are observed, including NiAs-type binary crystals, as well as the FCC 
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monocrystals of the large spheres. For sufficiently soft particles (Y < 50) the system remains 
in the fluid state at any pressure. Notice that Fig. [3] only reports the identity of the first 
crystal (the one formed at the lowest pressure) that nucleates from the fluid. In fact, as the 
system's Young's modulus is decreased and particles become penetrable a series of different 
crystalline structures can be generated upon increasing the external pressure. This peculiar 
behavior is intrinsic to soft particles [M] 125] , and the equilibrium phase diagram for spheres 
with this specific potential has been thoroughly described in ref. [26J. 

The pressure-structure dependence of these crystals is reported explicitly in Fig. [4] for 
s = 2,3 and 4, and can be found near the upper boundary of the diagrams where Y\ ~ Y 2 . 
Crystal structures were identified by combination of visual inspection, spherical harmonic- 
base order parameter q§ [35] , and by comparison with the known crystal structures reported 
in [3TH33] 136] [37]. The most abundant one is the NaCl-type, where the large spheres form 
an FCC lattice and the small ones fill the octahedral holes, as an interpenetrating FCC 
lattice. Next is the NiAs-type, where the large spheres pack in an HCP lattice and the small 
ones again fill the octahedral holes, this time in form of an interpenetrating simple hexagonal 
lattice. Coexistence of NaCl/NiAs-type crystalline phases in the AB-type mixtures has been 
observed experimentally and theoretically also in [321 [38] . The third crystal type is the FCC 
lattice of the large spheres, but no order for the small ones. The latter only develops for 
s > 3 and values of Y that are sufficiently small to allow hopping of the smaller components 
between different octahedral holes. Further reduction of the Young's modulus leads to the 
destabilization of any crystal structure as the system approaches the ideal gas limit. 

Figure [4] shows twelve two dimensional slices of the whole diagram, for a combinations 
of three size ratios and four different pressures. Depending on the combination of the four 
parameters we find a plethora of different phases, making the asymmetric case in which 
(Y 2 <C Yi) the more interesting one. In fact, a whole new class of phases becomes accessible 
to the asymmetric dumbbells, including spherical and cylindrical micelles, gyroid, lamellar, 
perforated lamellar and honeycomb phases. These phases are regularly observed in am- 
phiphilic and diblock copolymer systems [39] |4"T] . as well as in the polymer-nanoparticle 
tethered systems (401 H21 H2]. However, the aggregation in those systems is driven by both 
attractive and repulsive forces, while in our system no attractive force is present. 

The symmetry/nature of these mesophases can be understood by analogous arguments 
developed to describe the packing of amphiphilic molecules, i.e. the shape or more precisely 
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the effective shape of the dumbbells determines how they will pack in three dimensions. 
Indeed, shape asymmetry is key to the formation of mesophases. Shape asymmetry can be 
obtained in our system by either setting a different size for the two components constituting 
the dumbbell or by allowing one component to be softer than the other. Mesophases are 
observed when both physical size and softness asymmetry are present. We believe that this 
is because hard asymmetric dumbbells impose packing constraints that are too restrictive 
for them to arrange into curved surfaces without creating large stresses across the overall 
aggregate. This causes hard asymmetric dumbbells to form either crystalline or disordered 
phases. However, once the larger component of the dumbbell becomes deformable, packing 
constraints are relaxed and the formation of curved surfaces, with minimal interfacial energy 
and maximal entropy, becomes possible. Specifically, we find that, as Y 2 is decreased, the 
small and harder spheres can readily assemble into spherical clusters with overall FCC sym- 
metry or hexagonally-ordered cylindrical clusters; in both cases these clusters are stabilized 
and separated by a fluid layer formed by the overlapping larger spheres. 

As Y 2 is further decreased the large component becomes more deformable, effectively 
modulating the shape of the dumbbell, which becomes less conical. This leads to the forma- 
tion of first a gyroid phase of the smaller components, then a double gyroid phase (involving 
both components), and finally, as the effective shape of the dumbbell resembles more and 
more that of a cylinder, we find a lamellar phase. The thickness and the stacking of the 
lamellar phase depends on the specific values of Y%, Y 2 and pj. 

If we now also allow the smaller component to be deformed, additional curvature and 
degrees of freedom are introduced into the system which can now also form a perforated 
lamellar phase [44], obtained by interconnecting the soft layers through holes within the 
hard layers. 

Changing the size ratio of the dumbbell or modulating the relative strength of the Young's 
moduli of the two components is not the only route to phase change in these systems. 
As mentioned above, for Hertzian interactions, it is possible to induce the formation of 
different solid phases at constant s, Y\ and Y 2 by increasing the system pressure. Apart 
from the obvious transition from fluid to solid at low densities, we find that pressure can 
also drive several order-to-order transitions. For instance a gyroid phase can transform into 
a perforated-lamellar phase which can then turn into a lamellar phase and subsequently a 
honeycomb structure. Again, for large values of s the pressure does effectively alter the 
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overall geometric shape of the dumbbell by allowing for a greater degree of overlap among 
the larger components, which results into predictable symmetry changes of the aggregates, 
similarly to what we have previously discussed. 

The change of pressure introduces even more dramatic transformations in the Y\ ~ Y2 
region of the phase diagrams. Here we observe multiple re-entrant melting and crystal-crystal 
transitions upon pressure increase. Low pressure NaCl-type, NiAs-type and soft-sphere-FCC 
crystals melt under higher pressure and renucleate as CsCl-type binary crystals, or simple 
cubic (SC) and simple hexagonal (H) monocrystals of the soft spheres. At even higher 
pressures these crystals melt and again renucleate as A20-type monocrystals of the large 
spheres, where the small ones remain fluid, or HCP crystals of the large spheres, where the 
small spheres slightly overlap and aggregate into cylinders that fill the octahedral holes. 
Due to the bounded nature of the Hertz potential, many other crystal types would surely be 
found at even higher pressure, but we chose to stop compressing the system at p* = 1700. 

To test the robustness of our results, we also explored how structure formation is affected 
by the specific shape of the potential. Namely, we recompute the phase diagram for the 
case in which interparticle softness in described by a generalized inverse power law potential 
(GIPL). In this case, the extent of particle overlap is tuned by changing the exponent of the 
power law 

U{r lj ) = l LVW J (2) 

[0, Tij > Oij 

where = (<7j + <Jj)/2, and the exponent is defined as n,j = , which simplifies 
to the more conventional inverse power law form when rii = rij and U{ = <jj. In this case 
the hard-core limit is recovered for (n^ — > 00), and can be used to operationally define the 
diameter of the components to be equal to a^. Upon decreasing we can continuously 
allow for ever larger degrees of softness. Notice that, unlike the Hertz potential, the GIPL 
potential always presents an infinite barrier against complete particles overlap. 

For consistency with our previous results we restricted the computation of the phase 
diagram to <j\ < a 2 and n\ > 112- Specifically, we considered size ratios 2, 3 and 4, and 
for each s we explored several values of the exponents with n\ G [1,50] and 112 G [1,7]. 
This combination of exponents was selected in an attempt to find a set of parameters that 
are compatible to those of the Hertzian system in its most interesting region of the phase 
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diagram, i.e. we explicitly avoid the hard core limit for both components. All of our 
simulations are carried out using the same number of dumbbells at room temperature, and 
the system is equilibrated using the pressure-annealing scheme described earlier. 

Figure [5] shows the phase diagram of soft asymmetric dumbbells with GIPL potentials at 
three selected pressures and for the three different size ratios. Here is a general summary 
of what we observe: At low pressures all systems are in the fluid phase. Upon increasing 
the pressure, each system phase-separates into ordered aggregates. Similarly to what we 
observed for the Hertzian potential, the form of these aggregates depends both on size ratio 
and relative softness of two components. Our results reveal the existence of several ordered 
phases such as hexagonal cylindrical micelles, spherical micelles ordered into an FCC lattice, 
layered lamellar phases, perforated lamellar phase and a bicontinuous gyroid phase. Two 
types of the perforated lamellar phase are clearly distinguishable in different regions of the 
phase-space - those that with perforation through layers of the small spheres and those with 
perforation through layers of the large spheres. In Fig. [5] they are shown as separated 
phases. At intermediate pressures despite the overall system symmetry, both components 
can freely diffuse within their own phase, while as pressure is increased crystallization of 
both components takes place. 

Although the details of the phase diagram are different from those resulting from the 
Hertzian dumbbells, it should be noticed that using the GIPL potential we managed to 
reproduce almost all of the aggregates previously obtained. Due to our choice of exponents 
for the GIPL potential, a fair comparison between the two phase diagrams should only 
include relatively small values of Y2 in Fig. |4| The only exceptions are the crystal of the 
smaller components FCC1, and the nature of the gyroid phase (indicated between dotted 
lines in the GIPL diagram) which seems to be metastable. We checked this by running several 
simulations with non identical initial conditions but under the same set of thermodynamic 
parameters. We find that the gyroid phase does not consistently form as it instead does 
for Hertzian dumbbells. The metastability of the gyroid phase was discussed in detail for 
systems of Block Copolymers in jl5] and for polymer tethered nanospheres in |43j . A more 
specific study of this particular issue is out of the scope of our paper. Figure [6] shows 
snapshots of some of the crystalline phases found in this work, while Fig. [7] shows snapshots 
of the mesophases obtained. 

It is important to highlight that our results indicate that structure formation in these 
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systems, within the range of pressures considered here, is remarkably robust with respect to 
the specific choice of the potential, be that bounded or unbounded. Nevertheless, apart from 
the similarities in the phase structure, there is a fundamental difference between the two 
systems which is most evident when both components are soft. Namely, below a size-ratio- 
dependent value of Yi Hertzian dumbbells remain in the fluid phase no matter how large is the 
external pressure imposed on the system, while the fluid phase for the GIPL potential seems 
to always transforms into ordered aggregates once a sufficiently large pressure is achieved. 
For instance, we looked at the phase behavior of the GIPL system for s = 2 and n\ = n<i 
down to a value of the exponent of 0.3, and we indeed keep on finding formation of ordered 
aggregates. This is obviously not conclusive, as it is not clear whether there exists an onset 
value of the exponent n > (which could in principle be also size ratio dependent) below 
which the system stops crystallizing at any external pressure. This important but subtle 
issue is currently under investigation and the results will be published elsewhere. 

In conclusion, we have shown how the interplay between shape and deformability can 
be exploited to generate a large number of ordered mesophases in purely repulsive systems. 
We thoroughly explored the phase space for two cases of bounded and unbounded poten- 
tials covering a large range of systems spanning from the hardcore to the ideal gas limit 
(for Hertzian spheres). Our results show that the formation of most of the phases can be 
rationalized in terms of the effective geometry of the components and are not very sensitive 
to the specific choice of the potential. Some of the phases that we can easily obtain with 
soft nanoparticles, namely the gyroid, double gyroid and perforated lamellar phases, are 
also observed in block copolymers and have been very recently the subject of intense inves- 
tigation [331 Ho] . Nevertheless, bounded potential also give us access to multiple crystalline 
phases with pressure dependent symmetry, effectively providing for an almost unlimited va- 
riety of structures. Given the vast size of the phase space covered in this paper, we have not 
attempted to accomplish the heroic task of computing the equilibrium phase diagram and 
the relative stability of each of the many structures observed, but we hope that our results 
will inspire more theoretical and experimental work in this direction. 
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°2 

FIG. 1: Schematic representation of our model for a soft asymmetric dumbbell. Two soft spheres of 
diameter a\ and 02, with s = oijox > 1, are connected with a rigid bond of length d = 1/2(<ti+<T2). 
Both paricles are assumed to have equal mass. 
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FIG. 2: Illustration of the pair potentials used in this work to model the soft interactions between 
the dumbbellar components, a) The Hertz potential according to Eq. (1), plotted for the range 
of values of the Young's modulus 2 < Y < 10 4 . b) The generalized inverse power law potential 
(GIPL) according to Eq. (2), plotted for the range of the exponents 0.5 < n < 10. 
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FIG. 3: Crystallization of Hertzian dumbbells for size-ratios 1 < s < 3 and rigidities Y\ = Y% E 
[10, 10 6 ]. Below a size-ratio-dependent value of Y dumbbells stop crystalizing for any value of s 
(ideal gas limit). In the opposite regime (hard dumbbell limit) crystallization occurs below s = 1.05, 
where the spheres arrange into FCC lattice, and above s = 2.4, where the dumbbells arrange into 
an NaCl-type crystal. The width of the crystallization gap narrows as the repulsion become softer 
and other crystal types become accessible. Here we indicate the NiAS-type and crystals of the 
large spheres arranged into an FCC structure (FCC2), which are the ones that first nucleate from 
the fluid phase. Other structures can be formed upon further increase of the pressure (see text for 
details). 
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FIG. 4: Y\ — Y2 phase diagrams of Hertzian dumbbells for size ratios s = {2,3,4} and rescaled 
pressures p*j = {170,510,1020,1700}. Crystalline phases are obtained for Y% ~ Y% 3> 1, and 
include NaCl-type, NiAs-type and crystals of the large dumbbellar component, FCC2, which melt 
and recrystallize first into H and CsCl-type crystals and then into SC, A20 and HCP structures 
as the pressure is increased. For Y\ 3> I55 ~ 1 crystals of the small component, FCCl, are also 
observed. The region in between is dominated by mesophases including: FCC ordered spherical 
micelles (Sphere), hexagonally packed cylinders (Cy), gyroid and double gyroid (Gyroid), lamellar 
(Lamellar) and perforated lamellar (PL), honeycomb of the small components and cylinders of the 
large ones (HC1) and its inverted phase (HC2)15 
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FIG. 5: m — ri2 phase diagrams of soft asymmetric dumbbells with GIPL potentials, for size ratios 
s = {2,3,4} and rescaled pressures p*j = {86,136,1200}. Qualitatively, the same mesophases as 
in the case of the Hertz potential are obtained. The perforated lameral phase appears here as 
two distinguishable phases — perforated lamella through the layers of the large spheres (PL1) and 
perforated lamella through the layers of the small spheres (PL2). 
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FIG. 6: Some of the crystals obtained in the system of the asymmetric dumbbells. The large sphere 
is depicted in red (light color) and the small one is depicted in blue (dark color): (A) NaCl-type, 
(B) CsCl-type, (CI) and (C2) NiAs-type shown from two different perspectives, (Dl) 
and (D2) A20 crystal shown from two different perspectives. For the sake of clarity the 
sizes of the spheres has been altered. In snapshots (Dl) and (D2) the small spheres remain fluid 
and are not shown. 
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FIG. 7: Mesophases formed by a system of asymmetric dumbbells. Depicted in red (light color) 
is the larger and softer component of the dumbbell, while the smaller and harder one is depicted 
in blue (dark color). (A) FCC ordered spherical micelle -Sphere (B) hexagonally packed cylinders 
-HPC, (C) gyroid of the hard spheres, (Dl) and (D2) double gyroid phase -all referred to as Gyroid 
(E) perforated lamellar -PL, (F) lamellar, (G) honeycomb of the soft spheres and cylinders of the 
hard ones -HC1, (H) honeycomb of the hard spheres and cylinders of the soft ones -HC2. For the 
sake of clarity the larger spheres are portrayed with a smaller diameter. In addition, in snapshots 
(A), (B), (C) and (D2) the soft spheres are depicted using a light, low density pixel representation, 
while in snapshots (Dl) the same was done for the hard spheres. 
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